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Abstract 

We investigate the potential of the sparse data produced by the Catalina Sky Survey astrometric project (CSS for short) 
in asteroid shape and rotational state determination by the lightcurve inversion method. We show that although the 
photometric quality of the CSS data, compared to the dense data, is significantly worse, it is in principle possible that 
these data are for some asteroids with high lightcurve amplitudes sufficient for a unique shape determination. CSS 
data are available for ~ 180 asteroids for which shape models were previously derived from different photometric 
data sets. For 13 asteroids from this sample, we derive their unique shape models based only on CSS data, compare 
the two independent shape models together and discuss the reliability of models derived from only CSS data. We 
also use CSS data to determine shape models for asteroids with already known rotational period values, derive 12 
unique models and compare previously published periods with periods determined from the full 3D modeling by the 
lightcurve inversion method. Finally, we test different shape resolutions used in the lightcurve inversion method in 
order to find reliable asteroid models. 
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1. Introduction 

Orbital parameters are currently known for more than 
400000 asteroids. On the other hand, rotational states 
and shapes were determined only for a small fraction 
of them. In t he M inor Planet Lightcurve Databas^U 
( Warn er et all [2009), periods for ~ 3 500 asteroids 



are stored. Convex shapes with spin axis directions 
were derived only for ~ 200 asteroids, these three- 
dimensional models of asteroids are available in the 
Database of Aster oid Models f rom I nversion Tech- 
niques^ (DAMIT, iDurech et all 1201 Oh maintained by 
the Astronomical Institute of the Charles University in 
Prague, Czech Republic. 

All asteroid models in DAMIT were derived by the 
lightcurve inversion method (LI). This gradient-based 
method is a powerful tool that allows us to derive ba- 
sic physical properties of asteroids (the rotational state 
and the shape) from their disk-integrated photometry 
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(see 



2001 



Kaasalainen and Torppal 2001 : Kaasalainen et all 



2002). LI, in a form we use, leads to a convex 



shape approximation of a real asteroid (close to its con- 
vex hull) that has to be a single body (or the possible 
moon has to be so small that its photometric signature is 
indistinguishable from the noise) in a relaxed rotational 
state (i.e., rotates around the axis with a maximum mo- 
mentum of inertia). Sidereal rotational period and the 
direction of the spin axis are determined simultaneously 
with the shape. Typical number of free parameters is 
-50-80. 

We use two different types of disk-integrated photom- 
etry: (i) dense-in-time, which typically consists of tens 
to a few hundreds of individual data points observed 
during one revolution (typically several hours), and (ii) 
sparse-in-time, where the typical separation of individ- 
ual measurements is large compared to the rotation pe- 
riod. For sparse data, we usually have a few measure- 
ments per night, such as in the case of astrometric sky 
surveys. Based on the type of the photometry, we use 
the terms "dense lightcurves" and "sparse lightcurves". 
Unique asteroid shape models are typically feasible if 
we have dense lightcurves from at least three apparition, 
combined dense and sparse data give us more variabil- 
ity. 
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Our long-term strategy is to enlarge the number of 
asteroids with known shapes and rotational states, be- 
cause it can help us in understanding the physical pro- 
cesses that take place in the asteroid's populations, 
such as Near Earth Asteroids (NEAs), Main Belt As- 
teroids (MBAs) or even ast eroids in ind i vidual fami 
lies ( e.g., Koronis family, IShvaiu l2002t ISlivan et al. 



Response System, Hodappetall 20()3), and later also 

etaiTbOOll). and the 



from the Gaia satellite (IPerrvman et 



2003). Current distribution of periods and spin axes is 
the direct result of the evolution of these objects start- 
ing with their formation until the present time (several 
hundreds of Myr to ~ 4 Gyr for most studied aster- 
oids). The most prominent processes acting on aster- 
oids a re collision s and mass shedding, and the YOR P ef- 



fect dlubincaml Hjool |^^uhUckietafl|2003l). The 
knowledge of the shape can be used for several pur- 
poses, e;g^_JTjJntheconstruction of a thermal model 



(e.g., iMiiller et aU 1201 ll) . where values for geometric 



albedo, size and surface properties can be determined, 
(ii) a sample of real shapes instead of synthetic ones can 
be used for the statistical study of the non-gravitational 
forces ( YarkovskjQ and YORP effects), or (iii) in combi- 
nation with star occultations by asteroids. These events 
(observed for hundreds of asteroids) give us additional 
information about the shape (e.g. non-convexities) and 
can lead to a size estimation wit h a typical uncertainty 



of 10% (see Durechetal., 2011). 



Most of the currently available photometric data were 
already used in the LI. The only significant excep- 
tion are data from the Catalina Sky Survey astrometric 
project (CSS for short, Larson et ali 2003 ). In this pa- 
per, we investigate the possibility of determination of re- 
liable asteroid shape models from only a small amount 
of sparse-in-time distributed low-quality photometric 
measurements (~ 100) by the LI. Such data from CSS 
project are available for several thousands of asteroids. 
We want to find out if these data are for asteroids with 
high lightcurve amplitudes of a sufficient amount and 
quality for a unique shape determination, and if so, how 
reliable these asteroid models are. The investigation 
of sparse data capabilities and the reliability of derived 
models is important, because (i) it can lead to a deter- 
mination of new asteroid models without a need of ob- 
serving any additional photometric lightcurves, and (ii) 
in a few years, another huge amount of sparse data from 
three astrometric surveys will be available - from the 
Pan-STARRS (Panoramic Survey Telescope and Rapid 



3 Yarkovsky-0'Keefe-Radzievskii-Paddack effect, a torque 
caused by the recoil force from anisotropic thermal emission which 
can alter rotational periods and orientation of spin axes 

4 a thermal force acting on a rotating asteroid which can alter semi- 
major axes 



LSST (Large Synoptic Survey Telescope, llvezic et al 



2008). Understanding the CSS data with respect to the 
asteroid shape modeling will speed up future processing 
and use of the new sparse data. 

2. Photometry 

For a uniqu^l shape determination of an asteroid, 
we typically need dense photometric data from at least 
three apparitions and a good coverage of the solar phase 
ang 110 (i.e., different viewing geometries), ~ 20 such 
dense lightcurves are usually sufficient. The biggest 
limiting factor of these data is their insufficient amount. 
Dense lightcurves were observed for ~ 3 500 asteroids. 
Only in ~ 100 cases, asteroid shape models were de- 
rived. Determination of new models based on the dense 
data is now possible only if new lightcurves are ob- 
served. Many amateurs or semi-professionals use their 
telescopes for asteroid observations, but their main goal 
is the determination of the synodic period, which can 
be derived already from one apparition. If a period for 
a particular asteroid is then secure, they usually do not 
observe this object any more. We have to observe such 
object by ourselves. This approach is demanding on ob- 
servational time and also on fundings. Models based 
only on dense data are important in testing the relia- 
bility of models based on sparse data (due to possible 
comparison). 

Currently available sparse photometry is accessible 
via the Asteroids - Dynamic Site databas^l (AstDyS), 
where data from hundreds of astrometric observatories 
are stored. The photometry is mostly a by-product of 
astrometric measurements and in most cases, asteroid 
magnitudes are given to only one decimal place, so the 
accuracy is 0. 1 mag at best. Whether or nor this is suffi- 
cient for a unique shape determina tion for a rea s onabl e 
number of asteroids was studied in iHanus et al.l (1201 11) . 
The authors have found 7 observatories with quality 
data and used these data in combination with relative 
lightcurves for asteroid shape modeling. The most ac- 
curate sparse photometry is from the U.S. Naval Obser- 
vatory in Flagstaff (USNO-Flagstaff station, MPC code 



5 We define a unique solution as follows: (i) the best period has at 
least 10% lower \~ than all other periods in the scanned interval, (ii) 
for this period, there is only one pole solution with at least 10% lower 
X 1 than the others (with a possible ambiguity for A ± 1 80°), and (iii) 
this s olution fulfills all our additional tests (discussed in IHanus et all 
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689). These data were already studied in IDurech et al, 
(2009) and based on these data, the authors derived con- 
vex shape models for 24 asteroids. The biggest dis- 
advantage of the USNO-Flagstaff station data is that 
they are available only for about 1000 brightest aster- 
oids. Data from Hipparcos satellite jESAl 1 1997b and 
Roque de los Muchachos Observatory, La Palma (MPC 
code 950) are also accurate (compared to the other 
observatories), but their amount of measurements is 
even smaller than for the USNO-Flagstaff station. The 
largest amount of sparse photometry is available from 
the Catalina Sky S urvey observatory (CSS for short, 
Larson et al .1 l2003h. Th ese da ta were already used in 
the LI bv lHanus et al. ( 201 1 ). but only in combina- 
tion with other photometric data. Their typical photo- 
metric accuracy was ~ 12%. Although this accuracy 
of the CSS photometry seems very low, the data are 
valuable for asteroids with lightcurve amplitudes higher 
than > 0.3 mag, which is about a half of all aster- 
oids (based on the light curve data of the M inor Planet 
Lightcurve Database, IWarner et all 120091) . There are 
dense lightcurves for only ~ 700 asteroids with num- 
bers over ten thousand, on the other hand, we have at 
least 100 sparse data points from CSS for ~ 6 000 of 
these asteroids. 



3. Reliability of asteroid models derived only from 
the Catalina Sky Survey photometric data 

The Catalina Sky Survey astrometric project has pro- 
duced the largest amount of sparse photometric mea- 
surements. For most asteroids, no other data of suffi- 
cient quality are available. 

Test 1. One way how to test the reliability of mod- 
els based only on CSS data is to compare them with 
models derived from different photometric data set. 
So far, about a hundred of asteroid models were de- 



rived from dense data (e.g., 


Kaasalainen et all 


2002; 


To™ 


3a et al., 2003; Slivan et al., 2003; Marciniak et al., 


2007 


; Durechetal.. 20071). 


another ~ 100 models 



were derived from combined dense a nd sparse data 



(IDurech et al.L 12009; Han us et al.L 1201 II) . These models 
are believed to be reliable and thus correct. For most 
of these asteroids, CSS data are also available, so we 
tried to derive their models from only CSS data itself. 
We had at least 50 individual data points for 185 aster- 
oids. In 13 cases, we were able to derive their unique 
solutions. Only in 7 cases the solution was in an agree- 
ment with the model based on a different data set, other 
six solutions were formally correct, but derived periods 
were different from the expected ones. 



Test 2. The main inconvenient when computing mod- 
els from sparse data is that synodic periods are known 
only for ~ 3 500 asteroids (determined from dense 
lightcurves and stored in the Minor Planet Lightcurve 
Database). This a priori information is used in the shape 
modeling: we search for the model only near the as- 
sumed period (e.g., we test the initial period values from 
an interval (0.95 P, 1 .05 P), where P is the published 
period value). Typically, we test a huge number of ini- 
tial period values (usually thousands), because the mod- 
eled parameter space is full of local minims and the 
gradient-based LI method converges to the nearest local 
minimum corresponding to the initial period value. The 
minimal difference between two local minims is given 
mainly by the timespan of the observational data (see 
Eq. [1} and is usually very small (< 10" 3 h). Having a 
short period interval saves a lot of computational time 
and gives us confidence that the model is correct if the 
derived period value is close to P. Unfortunately, for 
most asteroids, we do not know the period (cannot be 
directly determined from sparse data without the model 
computation as in the case of the dense data), so we have 
to search for the model on a period interval of all possi- 
ble period values, typically 2-100 hours. This consumes 
a lot of computational time and we loose the possibility 
of comparing the two period values. 



In this test we used the a priori information about the 
periods: for some asteroids, we knew their periods, we 
had only CSS data (there exists their dense lightcurves, 
but we do not have them in an electronic form), and, 
based on that data, we were also able to derive their 
models near the published period. We extended the 
model search on the period interval of 2-100 hours and 
tried to find the model solutions again (as would be done 
if the period was unknown). This was performed for 
12 asteroids. In 7 cases, the previous solution derived 
on a shorter period interval was reproduced also on a 
larger period interval. For 4 asteroids, we did not get 
a unique period solution on the extended period inter- 
val, and a formally correct (i.e., it fulfilled our condi- 
tions for a unique solution), but a different solution was 
found for one asteroid. On Fig. [1] we show a peri- 
odogram of asteroid (5647) 1990 TZ, where for each 
initial period value a^ 2 -value corresponding to the best 
shape model and pole direction (a local minimum in 
the multi-dimensional parameter space) is plotted. This 
model computation was based on 87 individual mea- 
surements from CSS, derived period P = 6.13867 h is 
in agreement with a perio d P = 6.141 h reported by 
Bembric kaiid Bolt! d2003h . 
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Figure 1: Periodogram of asteroid (5647) 1990 TZ. Each point corresponds to a particular local minimum in the modeled parameter space of 
periods, pole directions and shapes. 



Test 3. Since we derived several unique models that 
were clearly incorrect, we performed this test to check 
more carefully the stability of the solution and to de- 
tect false solutions. The shape of an asteroid is repre- 
sented by coefficients of its expansion into the spherical 
harmonic functions. These coefficients are optimized in 
the LI. The number of coefficients is given by the order 
n of the shape expansion we use, we call n the shape 
resolution. In the LI, we typically use n — 6, which 
corresponds to 49 coefficients ((n + l) 2 ). For asteroids 
with < 100 sparse data points from CSS, so many co- 
efficients could be to much and the LI could then lead 
to a random unstable solution that is even unique and 
physically correct. Another input parameter of the LI 
is the period step P step between two subsequent periods 
from the scanned interval. The minimal difference be- 
tween two local minims in the modeled parameter space 
of periods is given by 



P 1 

AP = — , 
T 



(1) 



where P is the period value and T the timespan of the 
observational data. In our computations, we use for the 
period step a value P step = 0.8AP. 

We tested the stability of the solution on different val- 
ues of the resolution n and on two values of the pe- 
riod step P st£p (0.8AP and 0.5AP). Table Q] summa- 
rizes the results of this test for all asteroids in Test 1 
for which we derived their unique solutions for n — 6 
and P ste p = 0.8AP (corresponds to column 6). Correct 
unique solution is marked by wrong unique solu- 
tion by "x" and no unique solution by "— ". 



4. Discussion 

In Test 1, where we compared the derived models 
with corresponding models based on different data sets 
(stored in DAMIT), in 7 cases out of 13, the solution 
was in a well agreement with a model based on a differ- 
ent and higher data set. From this point of view, CSS 
data, which represents the largest dataset of sparse pho- 
tometry with photometric accuracy theoretically suf- 
ficient for a unique shape determination, seem to be 
promising for the asteroid shape modeling. Unfortu- 
nately, several clearly wrong models that fulfilled our 
conditions on a unique solution were derived (they had 
different sidereal rotational periods than was expected). 
Similar problem occurred in Test 2, where we compared 
the derived models with their previously known periods. 
One model was clearly incorrect. 

There are several possible explanations why we de- 
rived incorrect models from the CSS data: (i) data for 
asteroids with low amplitudes are more noisy, (ii) low 
amount of data with respect to the number of modeled 
parameters, (iii) systematic errors of the survey. In cases 
(i) and (ii), the data are poor, contain not enough infor- 
mation about the period and shape, and could produce a 
random solution with an unrealistic rms value (the fit 



is "too perfect"). lOszkiewicz et al. ( 201 1 ) show that 
photometric data sets from astrometric stations (such 
as CSS) exhibit magnitude variations with apparent V- 
band magnitude. For the brightest asteroids, the images 
are saturated, and for the faintest, the background sub- 
traction is imperfect. This consequence of (iii) is prob- 
ably partially visible in our data - most of the incorrect 
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Table 1: Test of the stability of the model solution on the shape resolution n and the step in the period P ste p- First column gives the asteroid number, 
symbol "/" means that f or a given n and Pgtepj 

both unique and correct solution was derived, symbol "x" means that a unique but wrong solution 
was derived, and finally, in the case if no unique solution was determined we use the symbol "— ". 
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models are brighter (have higher absolute magnitudes) 
than the correct ones. We are not able to distinguish be- 
tween the effects of the amplitude and the noise from 
each other only from the sparse data itself without a 
model computation, so there is no way how to detect 
the case (i). The amount of data sufficient for a cor- 
rect unique shape determination is dependent mainly on 
the time epochs (geometry of observations) and also on 
the signal-to-noise ratio, which we do not know. So, 
the sufficient amount of sparse measurements strongly 
varies from one asteroid to the other. 

As seen from Table [TJ models derived correctly for 
resolution n — 6 and period step P step = 0.8AP are 
mostly determined for other resolutions n and both pe- 
riod steps P st ep. No wrong solution appeared. In several 
cases, a unique solution was not found, but the correct 
period was still detected (we had multiple pole solu- 
tions). On the other hand, for incorrect models from 
Test 1, we got unique solutions for different resolutions 
n and period steps P step only rarely and never for n — 2 
or n = 3. All these unique solutions were always also 
incorrect. There is a weak dependence on the period 
step, P step = 0.5AP seems to be more convenient be- 
cause when used, less incorrect models were derived. 
To minimize the determination of incorrect models that 
satisfy our requirement on unique solutions, we should 
use the values n — 6 and P step = 0.5AP for resolution 
and period step, respectively. If a unique solution is then 
determined, the computation should be rerun with n = 3 



and P step = 0.5AP. The model is then reliable if we get 
similar unique solutions. 

We rerun the computation for all 8 successfully de- 
rived models from Test 2 with n = 3 and P step = 0.5AP. 
Previously correct solutions were confirmed, the incor- 
rect one was ruled out. 

Our method assumes a single body in a relaxed ro- 
tation state (rotates around the principal axis with the 
maximum momentum of inertia). Objects like tumblers 
or binary asteroids are not identified by the used form 
of the LI, their model determination usually fails be- 
cause the single period model does not cover the more 
complex period features. The only exception are fully 
synchronous close binary syste ms with similar-sized 
components, such as 91 Antiope (IMerline et al. , 2000). 
These objects are reproduced as highly elongated single 
bodies, and because they have typically high lightcurve 
amplitudes, they could be easily present in the possible 
sample of models based on CSS data. 



5. Conclusions 

Sparse data from the Catalina Sky Survey are low ef- 
ficient in unique model determination and seem more 
accurate for less brighter asteroids. Incorrectly deter- 
mined models could be indicated by unrealistic rms val- 
ues, but not generally. 

Correct period solutions for a single asteroid but dif- 
ferent resolution n were always stable, the dispersion in 
the period was typically lower than AP (given by Eq.[T]i, 
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which is for a period value of 10 hours and the timespan 
of the observations 10 years ~ 10~ 3 h (typical values 
for the period and the timespan of the CSS data). The 
dispersion in the ecliptic latitude of the pole direction 
was usually < 20°, but in the ecliptic longitude of the 
pole direction often > 50°. This significant dispersion 
in longitude was partially caused by models with high 
values of ecliptic latitudes (longitudes are more dense 
near the pole than near the equator) and because of the 
small resolution of the shapes. 

We showed that CSS data can be successfully used 
for reliable asteroid model determination. The most ac- 
curately derived parameter is the rotational period. Al- 
though the Gaia satellite will produce during its oper- 
ational time in average only ~ 60 data points for each 
asteroid (for CSS, we have ~ 100), the photometric ac- 
curacy will be incomparably better than for CSS data. 
Less data will be compensated by better accuracy al- 
lowing determination of many new asteroid models (or 
at least accurate values for rotational periods). The com- 
bination of Gaia data with data from Pan-STARRS and 
LSST projects will significantly increase the amount of 
new models. For asteroids with low amplitudes, their 
signal-to-noise ratio of the Gaia data will be also low 
and we will encounter the same problem of incorrect 
solutions as for CSS data. The detection of these false 
models will be similar. 

Finding an infallible test for detecting the false so- 
lutions derived from Catalina data could allow us de- 
termination of shapes and physical parameters for tens 
of asteroids. While the LI is more successful for as- 
teroids with higher lightcurve amplitudes, the derived 
sample of asteroid parameters will be significantly bi- 
ased. For example, shapes will be more elongated than 
it is usual in the whole population. Careful de-biasing 
of the model sample will be necessary before its use in 
any statistical study and physical interpretation. 

Possible models based on Catalina data are good can- 
didates for follow-up observations, new dense photom- 
etry should help in the detection of incorrect solutions 
and should lead to more detailed shape models. 
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